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DESIGNING FOR IMPACT RESISTANCE WITH UNIDIRECTIONAL 

FIBER COMPOSITES 

by Christos C. Chamis, Morgan P. Hanson, and Tito T. Serafini 

Lewis Research Center 

SUMMARY 

Composite micromechanics, macromechanics, and the miniature Izod impact test 
are used to investigate the impact resistance of unidirectional composites. Several 
composite systems are examined both theoretically and experimentally. The composites 
are classified theoretically relative to their impact resistance for longitudinal, trans- 
verse, and shear modes. Experimental results are reported only for Izod impact with 
the fibers either parallel or transverse to the cantilever longitudinal axis. Impact re- 
sistance design criteria which evolved during this investigation are used to design hybrid 
composites with improved impact resistance. This is illustrated theoretically and dem- 
onstrated experimentally. Approximate design procedures using the impact factor are 
described. 

The effect of microresidual stress on the longitudinal impact resistance is examined 
theoretically for composites with fiber-to-matrix modulus and/or stress ratio of approx- 
imately four. 

Predicted and measured results are compared on a rank or order basis. Measured 
results consist of those obtained in the experimental portion of this investigation and 
those available in the literature. Scanning electron microscope photomicrographs of 
specimen fracture surfaces are also included. Experimental results are presented to 
sh,ow the variation of the transverse impact resistance as a function of composite 
intralaminar- shear -strength. Photographs of the various impacted test specimens are 
presented to illustrate the types of failure. 

The results show that the in-situ fiber elongation-to-fracture controls longitudinal 
impact. Debonding and delamination are controlled by matrix modulus and in-situ ma- 
trix elongation-to-fracture. Microresidual stresses are detrimental in fiber/metal ma- 
trix composites. The ranking of predicted and measured impact resistance is in excel- 
lent agreement for several composites which had been evaluated by various methods. 

The experimental results indicate that impact can result in three main types of 
composite fracture modes. These are cleavage, cleavage with fiber pullout, and de- 
lamination. Combinations of these modes also take place. The impact resistance of one 
hybrid composite system investigated was greater than either of the constituent 
composites. 



INTRODUCTION 


An important design aspect of fiber composite structural components is their impact 
resistance. Some basic work on impact resistance and on other closely related proper- 
ties of these materials has been reported in the literature. See, for example, refer- 
ences 1 to 5. However, the understanding of impact resistance of fiber composites has 
not advanced to the point where components can be designed for impact using conven- 
tional design procedures. 

To obtain an insight into the impact resistance of structural components made from 
fiber composites, we begin by examining their physical make-up. The components con- 
sidered herein are made by laminating several plies; the ply is itself a unidirectional 
composite. A better understanding of component impact resistance can then be obtained 
by investigating the impact resistance of individual plies, multilayer unidirectional com- 
posites, the interply matrix layers, and the constituent material properties and fabrica- 
tion processing variables. This report deals with such an investigation. The investiga- 
tion is limited to gross-type impact (sufficiently long impact contact times so that the 
entire component resists the impacting force) and to unidirectional composites which ex- 
hibit a linear static stress-strain relation to fracture. 

The objectives of the investigation are to obtain a better understanding of impact 
resistance through elementary theoretical considerations and simple experiments. The 
experiments are of a qualitative nature and serve as a means to rank the composites. 

The following factors are examined: interpretation of impact resistance in terms of 
the energy under the static stress-strain diagram; relation of this energy to constituent 
material properties and fabrication processing variables; identification of prevalent 
failure modes; identification of constituent material properties which have a strong in- 
fluence on impact resistance; construction of design criteria for improving impact re- 
sistance; and classification of several available fiber composites on an impact resist- 
ance scale. 

The theoretical expressions for predicting impact resistance are covered in the 
section THEORETICAL INVESTIGATION. Here, impact resistance associated with 
single or combined fracture modes is presented and discussed. Design concepts using 
hybrid composites and the impact factor are also covered. The detailed derivations are 
given in the appendix. The experimental investigation is described in the section EX- 
PERIMENTAL INVESTIGATION. In this section, the constituent materials, fabrication 
process, test specimens, and test methods are described. The experimental results 
are also discussed in this section. Both theoretical and experimental results are pre- 
sented in tabular and graphical forms and can serve as an aid in design. 
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SYMBOLS 


A cross-sectional area 

Ap delaminated area 

a constant defined in eq. (5) 

B constant defined in eq. (5) 

b width 

C-^ correlation constant in eq. (8) 

dj fiber diameter 

E modulus 

G shear modulus 

g gravitational constant 

H height weight dropped 

h member depth 

IED impact energy density 

IF impact factor 

K spring constant 

Kjp volume ratio of pullout fibers 

Kj c fracture toughness opening mode 

K nc fracture toughness shear mode 

kj fiber volume ratio 

ky void ratio 

L c member length over which uniform stress exists 

l length 

l cr fiber debonded length 

delaminated length 
Njp number of pullout fibers 

N T n number of delaminated layers 

S unidirectional composite (ply) strength 

£>* modified S, eq. (17) and (A56) 
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SfT fiber strength 

S®nc longitudinal compressive strength 

T temperature 

AT temperature difference between composite processing and use temperatures 
U energy, strain energy 

V volume 

v impacting weight velocity 

W impacting weight 

x, y, z structural axes coordinate system 

1, 2, 3 material axes coordinate system 

a thermal coefficient of expansion 

/3 correlation coefficients 

)3 void strain magnification on in-situ matrix 

e strain 

e* composite limit fracture strain 

c p matrix strain-magnification-factor 

H' 

<y stress 

r shear strength for interface bond 

Subscripts: 

C compression 

c core 

cr critical 

D debonding, delamination 

FPO fiber pullout 

f fiber property 

i summation index 

L longitudinal 

l unidirectional composite (ply) property 

m matrix property 

mp matrix limiting property 
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R residual stress 

S shear 

s shell 

T tension 

x, y, z directions coinciding with structural axes 
1, 2, 3 directions coinciding with material axes 
Superscripts: 

a averaged properties 

c core composite 

s shell composite 


THEORETICAL INVESTIGATION 


In general advanced unidirectional fiber composites exhibit linear stress-strain be- 
havior (fig. 1). Linear stress- strain relations are also retained at high rates of load- 



Figure L -Typical stress -strain curves of unidirectional 
fiber composite material subjected to high rate of loading. 
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LONGITUDINAL TRANSVERSE SHEAR 

(A) UNIAXIAL IMPACT. 



LONGITUDINAL 

(B) CANTILEVER IMPACT. 


TRANSVERSE 


Figure 2. - Composite geometry and impact loadings. 


ings (ref. 6). These linear stress-strain relations and composite micromechanics 
(refs. 7 and 8) form the basis of the theoretical development for computing the impact 
resistance. 

The impact loadings, which are considered here, are illustrated in figure 2. As 
can be seen in this figure, the impact loadings are either along the material axis of the 
composite (longitudinal, transverse, or shear) or at the free end of a cantilever. 


Longitudinal Impact Resistance 

Longitudinal impact loading can result in either of two modes of fracture; these are 
(1) cleavage - the fracture surface consists of fractured fibers and matrix which lie ap- 
proximately in the same plane and (2) cleavage with fiber pullout - the fracture surface 
consists of fractured fibers in combination with debonding and fiber pullout. In the lat- 
ter case not all of the fracture surfaces of the fibers lie on the same plane. Both of 
these fracture modes are extensively discussed in references 3, 4, and 9. 


Impact -Induced Cleavage Fracture 

The equation describing cleavage failure due to impact is obtained by determining 
the strain energy density. It is shown experimentally in reference 10 that the strain 
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energy density correlates with Izod impact. For longitudinal impact (fig. 1(a)), this is 
simply 


or 



e mT S £llT V 


( 1 ) 



( 2 ) 


where U is the strain energy, e* is the fracture strain, S is the fracture strength, V 
is the volume, and E is the modulus. The subscript group Z11T is defined as follows: 
l refers to unidirectional properties, 11 identify outward normal to the plane and 
stress directions in that order, and T identifies the sense of the stress. By using 
composite micromechanics (ref. 8) two equations can be derived for Sj 11T depending 
on whether the fibers or the matrix offer the primary resistance to fracture. The de- 
rivations are given in the appendix. Only the final equations are given here. The im- 
pact energy density (IED) equals the strain energy divided by the volume. The IED of 
composites with an Ej/E m ratio greater than 20 is approximated by 


IED - (1 " k v )kf/3 fT S fT (; 

2E f 

with an approximation error of less than 5 percent. The undefined variables in equa- 
tion (3) are as follows: ky and denote void and fiber volume ratios, respectively; 
/3j,p represents the in- situ fiber strength efficiency which reflects the fabrication pro- 
cess. The subscript f refers to fiber property. The important points to be noted in 

equation (3) are the quadratic dependence of the strain energy density on the fiber 
2 2 

strength Sjrp and the fabrication process variable For a high impact resistance 

composite, equation (3) imposes the following requirements: a high strength low modu- 
lus fiber, approximately 100 percent fiber properties translation efficiency, high fiber 
volume ratio, and low void volume ratio. Three additional points to be noted here are 
the following: 

(1) The dependence of the strain energy density and therefore impact resistance on 
SfT/Ef and kj has been clearly demonstrated in references 11 and 12. 

(2) The contribution of (1 - kj)/3^p is contradictory to the results predicted by the 


7 



5 



I 

.06 


Figure 3. - Potential impact resistance of fiber composite materials from table I. 


TABLE I. - LONGITUDINAL IMPACT RESISTANCE OF VARIOUS 
FIBER/RESIN MATRIX COMPOSITES 


[Fiber volume ratio k^ = 0. 5; void ratio k^ = 0; in- situ fiber 
strength efficiency ( 3 ^ = 1. 0.] 


Fiber 

Density 

j Fiber modulus 

Fiber strength 

Predicted longitudinal impact 


, 3 

„ 3 

2 


_ 2 


| energy density 



gz cm 

lb/ in. 

N/cm 

psi 

N 7 cm 

ksi 











cm -N/cm ^ 

in. -lb > in. ** 

Rank 

Boron 

2. 62 

0. 095 

414X10 5 

60X10 6 

3. 18X10 5 

460 

608 

880 

7 

E -glass 

2. 49 

. 090 

69 

10 

2. 58 

360 

2240 

3250 

3 

S-glass 

2. 49 

.090 

86 

12. 4 

4. 62 

670 

6250 

9050 

1 

Modmor-I 

1. 99 

. 072 

414 

60 

1. 73 

250 

179 

260 

10 

Modmor-II 

1. 74 

. 063 

262 

38 

2. 52 

350 

656 

950 

6 

ThorneI-400 

1. 71 

.062 

207 

30 

2. 90 

420 

1020 

1470 

5 

Thornel-50 

1.63 

.059 

345 

50 

1. 66 

240 

200 

290 

9 

Thornel-75 

1. 85 

.067 

473 

75 

2. 62 

380 

324 

470 

8 

PRD-49 

1. 38 

.050 

173 

25 

2. 76 

400 

1040 

1600 

4 

UARL-344 glass 

3. 60 

. 130 

128 

18. 6 

4. 83 

700 

4860 

7050 

2 
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debonding and fiber pullout mechanism. See section Longitudinal Impact with Cleavage 
and Fiber Pullout and also references 3 and 4. 

(3) Equation (3) is a simple and convenient means to rank fiber composites for lon- 
gitudinal impact resistance. 

A graphical representation of equation (3) for various available composites is shown 
in figure 3, where the strain energy density is plotted as a function of S^/E^ (ratio of 
fiber strength to fiber modulus) which equals in- situ fiber elongation-to-fracture. These 
same composites have been ranked according to equation (3) in table I. Note that in ta- 
ble I three relatively new fibers have been listed; they are Thornel-400, a high modulus 
organic fiber (PRD- 49), and UARL-344 glass (ref. 13). 

Rank comparisons of results reported in the literature with those predicted by 
equation (3) are shown in table n for notched Charpy impact, in table in for fracture 


TABLE n. - COMPARISON OF PREDICTED RESULTS WITH NOTCHED CHARPY IMPACT DATA (REF. 13) 


WITH FIBERS PARALLEL TO LONGITUDINAL AXIS OF BEAM 


Composite 

Fiber 

volume 

ratio 

Fiber modulus 

Fiber strength 

Impact energy 

N- cm 2 

psi 

N cm 2 

ksi 

Measured (ref. 13) 

Predicted a 

cm-N 

ft -lb 

Rank 

cm-N' cm 

in. -lb in. 

Rank 

Thornel-50 epoxy 

0. 55 

345:<10 5 

50xl0 6 

1. 66xlO S 

240 

544 

4 

4 

218 

315 

4 

Boron epoxy 

. 55 

404 

58. 5 

3. 18 

460 

1356 

10 

3 

687 

995 

3 

UARL-344 glass epoxy 

. 633 

128 

18.6 

4. 83 

700 

4080 

30 

2 

5760 

8 350 

2 

S-glass epoxy 

. 65 

85. 5 

12. 4 

4. 62 

670 

7340 

54 

1 

8160 

11 800 

1 


a Void ratio k v = 0; in-situ fiber strength efficiency = 1. 0. 


TABLE III. - COMPARISON OF IMPACT ENERGY DENSITY WITH MEASURED 


FRACTURE TOUGHNESS FOR GLASS-FABRIC COMPOSITES 


Specimen a 

Fiber 

Measured fracture toughness. 

K b 
K Ic 

Predicted impact energy 


volume 

N/cm 2 -cm 1 ^ 2 

ksi-in. 1 / 2a 


density 



ratio a 

Rank 

cm-N/ cm^ 

in. -lb/ in. 2 

Rank 

1 - 18 

0. 545 

27. 2 to 28.8 

24.8 to 26. 2 

3 

1700 

2460 

3 

2 - 18 

. 476 

22. 4 to 25.0 

20. 4 to 22. 7 

4 

1460 

2120 

4 

3 - 18 

. 589 

40. 8 to 41. 8 

37. 2 to 38. 1 

2 

1830 

2650 

2 

00 

T— ( 

1 

. 676 

43. 1 to 46. 8 

39. 2 to 41. 7 

1 

2110 

3050 

1 

5-18 

. 294 

17.4 to 20.6 

15.8 to 18.7 

5 

912 

1320 

5 

6-18 

. 225 

16. 4 to 18.0 

14. 9 to 16. 4 

6 

701 

1015 

6 


a Results from table 6 (ref. 11). 

^Results were obtained from beam splitting tests. 

c Fiber strength S fT = 207 N/cm 2 (300 ksi); fiber modulus E f = 6. 9xl0 6 N/cm 2 (10xl0 6 psi); 


void ratio = 0; and in-situ fiber strength efficiency = 1. 0. 
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TABLE IV. - COMPARISON OF IMPACT ENERGY DENSITY WITH MEASURED DATA OF GRAPHITE 


FIBER/EPOXY COMPOSITES AT CRYOGENIC TEMPERATURES 


[Temperature, T = 20. 3° R (-423° F), ref. 19.] 


Matrix type (measured) a 

Fiber 

Fiber modulus 13 

Fiber strength** 


Composite fracture energy 



2 

cm-N/cm 

9 

in. -Ib/in . 1 

Rank 

volume 

N/cm 2 

psi 

N/cm 2 

ksi 

Measured 3- 


Predicted impact energy density 









2 

cm-N/cm 

in. -lb/in. 2 

Rank 

q 

cm-N/cm 

in. -lb/in. 2 

Rank 

7. 3 

10.6 

3 

0. 38 

324X10 5 

47X10 6 

1. 59X10 5 

230 

46 

67 

4 

148 

214 

4 

7. 1 

10. 4 

4 

.39 

338 

49 

1.64 

238 

29 

42 

5 

155 

225 

3 

2.9 

4. 2 

6 

.40 

345 

50 

2. 54 

368 

152 

220 

1 

372 

540 

1 

8.2 

11. 9 

2 

. 40 

297 

43 

2.68 

289 

64 

93 

2 

279 

390 

2 

8.5 

12. 3 

1 

.42 

332 

48 

1. 50 

218 

58 

84 

3 

144 

208 

5 

3. 4 

5.0 

5 

. 41 

311 

45 

.99 

143 

L 23 . 

33 

6 

64 

1 93 

6 


a Measured by beam splitting method (ref. 20). 
^Reported in ref. 20. 


toughness, and in table IV for energy-absorbed-to-failure at cryogenic temperatures. 
As can be seen in these tables, the ranking comparisons are in excellent agreement. 


Effects of Microresidual Stresses on Impact Resistance 

The contribution of the matrix to impact resistance is not negligible in composites 
having a strong and stiff matrix and a good interface bond. For these types of com- 
posites, usually Ej/E m < 10, which is typical for fiber/metal matrix composites. 

The governing equation for the impact energy density for this case is given by 

(4) 


( 5 ) 


The subscripts l , f, and m denote ply, fiber, and matrix properties, respectively; 
a is the thermal coefficient of expansion; and AT is the difference between the com- 
posite processing and use temperatures. 


IED = aE m IS mT - B 


-l 11 / 


where 


a = 


2E^ 


m 


B = AT((j f - o- m )E f E 


m 
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MATRIX CONTROLLED 

FIBER CONTROLLED 


20xl0 2 14 xl0 2 



V 


>- 
i s\ 


>- 

o 

cn 


o 

c 


16 


121 — 


4 ' — 


o 

? 

S 

CJ 

fc 



2 


0 .1 .2 .3 .4 .5 .6 .7 

FIBER VOLUME RATIO, k f 

Figure 4. - Theoretical longitudinal impact resistance of boron-silicon carbide/ 
titanium unidirectional composite. Processing temperature, 832° C (1500° F). 


One very important point to be noted in equations (4) and (5) is that the strain energy 
density depends significantly on the microresidual stress. The microresidual stress is 
represented by the parameter B in equation (4). This dependence has not been reported 
previously in the literature. It is suspected that the presence of microresidual stress 
in the matrix produced some of the trends reported in references 3 and 14. However, 
the authors of these references did not attribute the decrease in fracture energy to this 
phenomenon. 

The dependence of the strain energy density and therefore the impact resistance on 
the microresidual stress is illustrated in figure 4 for a boron-silicon carbide coated/ 
titanium system. Two sets of curves are plotted in this figure. One set is for matrix- 
controlled failure with and without residual stress. The other set is for fiber- controlled 
failure with and without residual stress. This second set was obtained from equations 
(4) and (5) by interchanging the subscripts f and m (see appendix). 

The important point to be noted in figure 4 is that impact resistance, or fracture 
toughness, is very sensitive to the presence of microresidual stresses. Therefore, 
interpretation of experimental results from composites with Ef/E m < 10 must take the 
microresidual stress into account. 

Longitudinal impact loadings resulting in partial cleavage failure with debonding and 
fiber pullout is a combined fracture mode. Description of this type of mode will follow 
the description of the single modes. 
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Transverse Impact Resistance 


Transverse impact loadings of unidirectional composites (fig. 2(a)) result in brittle 
fractures. The amount of energy absorbed to fracture during transverse impact is re- 
ferred to as the transverse impact resistance. The strain energy divided by the volume 
of the material is referred to as the IED. This IED as measured under the transverse 
stress-strain curve is shown in figure 1(b). The governing equation is derived from the 
stress-strain diagram in figure 1(b) and the micromechanics relations of reference 8. 
The detailed derivations are presented in the appendix. The resulting equation for the 
transverse IED is given by 


IED 


1 

2 



e mpT 

VV22, 


22 


( 6 ) 


The variables in equation (6) are the following: &22T * s corre l a ti° n coefficient re- 
flecting the fabrication process; e m p,p is the maximum transverse strain that the in- 
situ matrix will experience when the composite is loaded in the transverse direction; 

is the void magnification of the transverse matrix strain; cp ^22 * s matrlx 
transverse-strain-magnification factor which is a complex function of constituent moduli 
and fiber content; and E^2 is the transverse composite modulus. 

There are several important points to be observed in equation (6); they are the 
following: 

(1) The transverse impact resistance is a complex function of the fabrication pro- 
cess, material properties, and composite properties. 

(2) The degree of bond at the interface is reflected by ^22T’ the poorer the inter- 
face bond, the smaller the value for this coefficient. 

(3) Increases in either void or fiber content or both have inverse square effects on 
the transverse impact resistance. These effects result in more brittle composite be- 
havior. 

(4) The impact resistance increases linearly with the ply transverse modulus. 

(5) The impact resistance increases as the square of the in-situ matrix-fracture- 
strain. 

It is important to note that the in-situ matrix-fracture-strain is not the failure 
strain of the bulk matrix material. For nonmetallic matrixes the former is a small 
fraction of the latter (ref. 8). The difference between in-situ and bulk matrix-fracture- 
strain is not widely recognized. As a result, efforts to correlate theory with experi- 
ment and to develop matrix materials which would result in improved composite proper- 
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ties have usually failed. However, both of these disparities can be remedied with suit- 
able micromechanics models and appropriate experiments (results of current unpub- 
lished research by the authors). 

The graphical representation of equation (6) for typical fiber composites is shown 
in figure 5. In this figure the transverse IED has been plotted as a function of fiber 
volume ratio. Three important points to be noted in figure 5 are the following: 

(1) The impact resistance of graphite fiber/epoxy composites is insensitive to fiber 
volume ratio. 

(2) However, boron and glass fiber/epoxy composites become quite brittle at high 
fiber volume ratios (>0.65). 

(3) All fiber /nonmetallic composites have approximately the same impact resistance 
at about 0. 50 fiber volume ratio. 

The variation of the transverse IED as a function of matrix modulus, is shown in 
table V for a Modmor-I/epoxy composite. As can be seen in this table, the IED in- 
creases very rapidly with increasing matrix modulus. The two reasons for this rapid 
increase are (1) the matrix- strain-magnification factor <P ^22 decreases rapidly while 
the composite transverse modulus E^ 22 increases (table V) and (2) the fiber is aniso- 
tropic; that is, the transverse fiber modulus is about 0. 69x10® to 1. 3X10® newtons per 
square centimeter (1x10 to 2x10 psi). 
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TABLE V. - EFFECTS OF MATRIX MODULUS ON TRANSVERSE AN SHEAR IMPACT ENERGY DENSITY 


[Modm or -I/epoxy composite; fiber volume ratio = 0. 5; void ratio = 0.] 


Matrix modulus 

Strain magnifi- 

Unidirectional composite moduli 3 

Predicted impact energy dens 



cation factors 3 








N/cm^ 

psi 



Transverse 

Shear 

Transverse 

She 



Trans- 

Shear 










verse 


N/cni 

psi 

N/cni 

psi 

cm-N/cm 

in. -lb/in. 3 

cm-N/cm 3 

0. 69X10 5 

0. 1X10 6 

3.26 

4. 37 

1. 86X10 5 

0. 27X10 6 

0. 97X10 5 

0. 14x10® 

0. 93 

1.34 

4.3 

2.07 

.3 

1.96 

3. 54 

4.35 

.63 

2. 56 

.37 

5.67 

8.20 

17.0 

3. 45 

.5 

1.42 

2. 97 

5.86 

.85 

3. 80 

. 55 

14.6 

21. 10 

36.6 

4. 82 

.7 

1. 12 

2. 56 

6.89 

1.00 

5. 18 

. 75 

27. 5 

39.8 

61.9 

6. 89 

1.0 

1.00 

2. 12 

7. 95 

1. 15 

6. 14 

.89 

39. 7 

57. 5 

115 

10. 03 

1.5 

1. 00 ! 

1. 65 

9.05 

1. 31 

7.66 

1. 11 

45.2 

65. 5 

237 


in. -lb/in. 3 | 

6.2 

24.6 
53. 1 

89.6 
166 
343 


a Generated with computer code of ref. 18. 


Shear Impact Resistance 


Shear impact loadings of unidirectional composite (fig. 2(a)) result in relatively 
brittle fracture. The amount of energy absorbed to fracture during shear impact is 
called herein shear impact resistance. The corresponding IED as measured under the 
shear stress-strain curve is shown in figure 1(c). The detailed derivations of the gov- 
erning equation are described in the appendix. The resulting equation for shear is given 
by 


IED=l/ ' - 12sem P s ) 2 G. 12 (7) 

Note the similarity of equations (7) and (6). Corresponding terms have analogous mean- 
ings, namely: is corre l a ti° n factor; c m pg is the in-situ matrix shear- 

fracture-strain; /3 y is the void contribution to the matrix shear strain; <P 2 * s *be 
matrix shear- strain-magnification factor; and ^ is the composite shear modulus 
in the plane containing the fibers. 

The important points noted in discussing equation (6) apply to corresponding terms 
in equation (7) as well. One additional point to be noted is that equation (7) describes 
also intralaminar shear delamination which will be described subsequently. 

The graphical representation for typical fiber composites is shown in figure 6. In 
this figure the IED for shear is plotted as a function of the fiber volume ratio. The 
important points in figure 6 are the following: 

(1) Boron/epoxy composites are superior in shear impact as compared with other 
fiber/epoxy composites when the fiber volume ratio is less than about 0. 6. 
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Figure 6. -Theoretical shear impact resistance of unidirectional composites. 


(2) The shear impact resistance of isotropic boron and S-glass fiber/epoxy com- 
posites is very sensitive to fiber volume ratio. 

The variation of the shear IED as a function of matrix modulus for a graphite 
Modmor-I fiber/epoxy composite is shown in table V. As can be seen in this table, the 
shear IED increases very rapidly with increasing matrix modulus. The reason for this 
very rapid increase is the variation of the matrix shear -strain-magnification factor and 
the composite shear modulus (eq. (7) and table V) with increasing matrix modulus. It 
should be noted that the shear IED increases more rapidly than the transverse IED 
as can be seen by comparing corresponding columns in table V. 

It is interesting to compare equation (7) with equation (17) (from ref. 15) which was 
derived from fracture mechanics considerations. Equation (17) from reference 15 is re- 
peated here using the notation of this report for convenience 

K nc = C 1 S J12S 

where K nc is the critical stress intensity factor (fracture toughness) in the shear mode, 
is a correlation coefficient, S^g is the intralaminar (horizontal) shear strength, 
and ky is the void volume ratio. By using the micromechanics definition for -^S’ 
equation (17) from reference 15 can be expressed as 
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K 


nc 



Note that the fracture toughness in the shear mode Kjj c depends linearly on the in-situ 
matrix shear-fracture strain e mpS and the composite shear modulus. It is inversely 
proportional to the matrix shear-strain-magnification factor. The parameter 
Gj 12A 0 /il2 i s a nonlinear increasing function of the matrix modulus (table V). The 
matrix modulus then is a more important parameter in increasing Kjj c than the bulk 
matrix elongation-to-fracture. This observation is not widely recognized in the fiber 
composite research community. Note also that the shear impact resistance is more 
sensitive to the parameter e m ^g/ ( P ^2 than is the fracture toughness. The former de- 
pends quadratically on this parameter, while the latter only linearly. 


Longitudinal Impact Resistance from Fiber Pullout 

Fiber composite fractured surfaces usually exhibit some debonding and fiber pull- 
out. This fracture mechanism has been investigated extensively (refs. 1 to 4). 

The following two assumptions are made to derive the governing equation: (1) the 
energy absorbed during impact is expanded in pulling out the fibers and (2) the interface 
bond strength is approximated by the intralaminar shear strength. Assumption (2) was 
first introduced in reference 2. The detailed derivations leading to the governing equa- 
tion are given in the appendix. The result for the impact energy density from fiber 
pullout is given by 


IED = - (1 - k )(- AV 2 (9; 

8 \^12S e mpS/ G £ 12 

The symbols in equation (9) have been defined previously. Equation (9) describes IED 
due to fiber pullout as a complex function depending on fabrication process, fiber and 
void contents, constituent strength properties, and composite shear modulus. The 
variation of IED as a function of constituent elastic properties is not easily seen in 
equation (9) because the parameter <P ^yi^l 12 depends on fiber and void contents, and 
on the constituent properties in a complex way. This parameter is defined herein as 
the "debonding parameter” because it is an indication of the local interface shear bond. 
Its dependence on matrix modulus and fiber volume ratio is shown in figure 7 for 
Modmor-I fiber/epoxy composites. 
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Figure 7. - Debonding parameter for Mod mo r -I/epoxy composite with zero 
voids. 


The important points to be noted from equations (9) and (A9) in conjunction with fig- 
ure 7 are the following: 

(1) Local bonding is enhanced with increasing fiber volume ratios (up to about 0. 65) 
or increasing matrix modulus. 

(2) IED due to debonding can be increased by any or all combinations of the follow- 
ing: poor interface bond, low in- situ matrix elongation-to-failure, large 

ratio, and constituents selection which results in low composite shear modulus ^ 

It is important to note that the parameters which enhance IED from debonding and 
fiber pullout are quite detrimental to composite structural integrity with respect to 
static strength and stiffness. 


Impact Resistance Due to Delamination 

Delamination in the context used here refers to the delamination due to shear of 
interply layers in multilayered composites. The energy expanded is referred to herein 
as the "impact resistance due to delamination. " 

The governing equation to describe this resistance is based on the following assump- 
tions: 

(1) Delamination occurs when the interlaminar shear strength has been exceeded. 

(2) Several interply layers could delaminate simultaneously. The detailed deriva- 
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tions are given in the appendix. The resulting equation for the IED from delamination 
is given by 


IED = i G (10) 

2 \V#xl2 / 

where N LD is the number of delaminated interply layers. 

Note that equation (10) is identical with equation (7) except for the coefficient N T jy 
Therefore, the discussion following equation (7) and the important points noted there 
apply to equation (10) as well. The additional point to be noted from equation (10) is, 
that for improved impact resistance, design the part to assure multi-interply delamina- 
tion. This should be applicable to high velocity impact as well as low. 


Longitudinal Impact with Cleavage and Fiber Pullout 

This type of impact resistance results in fractured surfaces consisting of broken 
fibers with debonding and fiber pullout. It was referred to as cleavage with debonding. 
The governing equation is a combination of equations (3) and (9). The result for the 
impact energy density for this case is given by 


IED = (1 - k^) JI 
^ 2E f 


„2 , d f k fD l 

PfT K f 


v%12 


4L c \ /3 12S € mpS/ G i 12. 


(ID 


where L c is the length of the component subjected to uniform stress which causes fiber 
fracture. 

It is important to note that the fiber pullout contribution (second term in eq. (11)) to 
impact resistance in equation (11) is strongly dependent on L c> The following example 
will illustrate the point. Using typical values for a Modmor-I fiber/matrix composite 
and assuming 40 percent fiber pullout result in a contribution of approximately 0. 3/L c . 
This contribution is negligible for longitudinal impact where L c is quite large. How- 
ever, the fiber pullout contribution will be significant in the case of localized or bending 
impact. 
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The fiber pullout contribution will, in general, be negligible (less than about 1 to 2 
percent) if 



(12) 


Equation (11) indicates that composites with high fiber modulus and low intralaminar 
shear strength are good candidates for high impact resistance. Since equation (11) is a 
combination of equations (3) and (9), the discussion following these equations applies to 
equation (11) as well. 


HYBRID COMPOSITES TAILOR-MADE FOR IMPROVED IMPACT RESISTANCE 


Hybrid composite is the term used for a composite which consists of two or more 
different fiber/matrix combinations. Typical examples are Modmor-l/epoxy-glass/ 
epoxy-Modmor-l/epoxy, HTS/epoxy-Thornel-50/epoxy-HTS/epoxy, and others. 

Using these composites for improved impact resistance is a major contribution of 
this investigation. The concept was discovered during the experimental portion of the 
investigation. It was observed that some of the impacted cantilever specimens 
(fig. 2(b) longitudinal) exhibited combined fracture modes consisting of fiber breakage, 
fiber pullout, and interply delamination. 

The hybrid composite takes advantage of two or more of these modes to improve 
impact resistance. It is an important and useful concept in designing structural com- 
ponents in general. The impact resistance of hybrid composites is thus not a material 
characteristic. 

The concept is illustrated here, by applying it to the cantilever structure shown in 
figure 8. The governing equation for impact energy density is given by 


oS 

1 11T 
IED = - - 1X1 

2 F s 
E Z11 



16 N LD'^ s £ 12 


+ 7r / N fD d fV E fn\/ S fT \ ^ 

16 \ bh l A Q s /( Q s 
x / \ b Z12/ \ b ZllT/ 


where the superscripts a, s, and c represent averaged core-shell, shell, and core, 
respectively; the subscript 1 refers to unidirectional composite properties along the 
direction indicated by the numerical subscripts following Z; the variables b, h, and Z 
represent width, depth, and length of the cantilever, respectively (see fig. 8 also); d^ 
is the fiber diameter, is the number of fibers that pulled out, and N TF) is the 

number of layers that delaminated. 
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Figure 8. - Cantilever subjected to impact. 


Examining equation (13) reveals that the shear contribution depends on E ?ll/ G U2 

and that both fiber pullout and delamination depend on the parameter E ?ll/ G U2* This 
means that, in order to take advantage of the high shear contribution of fiber pullout 
and/or delamination, high longitudinal modulus, low shear modulus, and low intralam- 
inar strength composites should be selected. Some composites which meet this criterion 
are Thornel-50, Modmor-I, and PRD-49 fibers in a resin matrix. 

There are three other sets of parameters in equation (13) which need careful exam- 
ination in designing hybrid composites for improved impact resistance. These are 
(1) (h/Z)^ for the shear contribution, (2) h/N T j/ for delamination, and (3) dj? N^/bh^ 
for fiber pullout. 

The shear contribution will be greater than 3 percent when 



> 10 


E 


a 

11 1 



The contribution of the fiber pullout will be greater than 3 percent when 


> 0.02 

bhZ „s 

E l 11 


(A3 2) 


(14) 


The contribution of the delamination will be greater than 3 percent when 



( 15 ) 
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The following expression must be satisfied for delamination: 


12S 


•i© 


— mm 




UlT’ S ?llc) 


( 16 ) 


where the variable denotes longitudinal compressive strength. 

Equation (13) in conjunction with the inequalities (eqs. (14) to (16) and (A32)) pro- 
vides relations which can be used to select parameters in designing composites with im- 
proved impact resistance. These equations were used in this investigation to guide the 
selection of the hybrid composites. 

The inequalities (eqs. (A32) and (14) to (16)) can be expressed in terms of constitu- 
ent properties by using the micromechanics relations for j^g, an d 

E Z11- 


APPROXIMATE DESIGN FOR IMPACT 

Structural components subjected to impact are designed using an equivalent static 
load. This type of design is the strength of materials approach and is a first-order ap- 
proximation (ref. 16, ch. 5). 

The governing equations are 


where 



(17) 



(18) 


and where is the allowable stress to account for impact, IF is the impact factor, 
is the static composite strength, K is the spring constant which depends on the type of 
impact (Relations for K for the types of impact described previously are given in ta- 
ble VI.), W is the impacting weight, H is the height from which W is dropped, and 
v /2g is the potential energy of W. 

Note the subscripts in and depend on the type of impact loading. For lon- 
gitudinal impact, for example, the subscripts will be HT(S^jjrp and ^u T ). The 
corresponding spring constant is given in the first line of table VI. 


TABLE VI. - SPRING CONSTANTS RELATIONS FOR 


VARIOUS IMPACT LOADINGS 


Type of impact loading 


Spring constant, 


Longitudinal 


Transverse 


Shear 


K 



Cantilever longitudinal (rec- 
tangular section including 
shear contributions) 


E Ul bh 


4 1 ' 


1 + 

Pi 12 


Cantilever transverse (rec- 
tangular section including 
shear contributions) 


E i 22 bh 


1 + 1. 2/-)Y^t22 

G / 23 


The moduli appearing in the spring constant relations in table VI can be evaluated 
using micromechanics. Values of these moduli as functions of void and fiber volume 
ratios are given in reference 7 for several composites. The computer code described 
in reference 17 may be used to generate additional ones. 


EXPERIMENTAL INVESTIGATION 

This portion of the investigation consisted of carrying out miniature Izod (ref. 18) 
impact tests to verify qualitatively the theoretical considerations and concepts described * 

in the THEORETICAL INVESTIGATION. 

i 

Materials and Specimen Fabrication 

Graphite, glass, and PRD -49 fibers in an epoxy resin matrix were used in the exper- 
imental investigation. The various fibers are listed in table VII. All fiber material 
was drum wound and impregnated with the epoxy resin ERL 2256-ZZL0820 (27. 0 pph 
resin). 
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TABLE Vn. - MINIATURE IZOD IMPACT DATA FOR FIBER/EPOXY a COMPOSITES 


Fiber 

Type 

Surface treatment 

Fiber 

Average impact energy 


Rank 





volume 












ratio 

Longitudinal 

Transverse 

Longitudinal 

Transverse 





cm-N 

in. -lb 

cm-N 

in. -lb 

Meas- 

Predic- 

Meas- 

Predic- 









ured 

ted 

ured 

ted 

Graphite 

Thornel-50S 

(b) 

0. 532 

85. 9 

7.6 

7. 9 

0.7 

5 

5 

3 

3 


Thornel-50 

Polyvinyl alcohol 

. 583 

208,0 

18.4 

3. 4 

. 3 

4 

4 

5 

5 


HTS 

(c) 

. 523 

56. 5 

5.0 

14. 7 

1.3 

6 

6 

2 

2 


Modmor-I 

None 

. 542 

215.0 

19.0 

4. 5 

.4 

3 

3 

4 

4 

Glass 

S 

(d) 

0. 486 

757. 0 

67.0 

15. 8 

1.4 

1 

l 

1 

1 

PRD-49 


— - 


280.0 

24. 8 

3. 4 

0. 3 

2 

2 

5 

-- 

Graphite 

HTS/Thornel-50S 

— 

0.598 

116. 3 

10. 3 

11. 3 

1.0 

3 

— 

— 

— 


HMS/Modmor-I 

--- 

.536 

132.0 

11. 7 

— 

--- 

2 

— 

— 

___ 


HMS 

— 


232.0 

20. 5 

— 

--- 

1 

--- 

--- 



a Epoxy resin: ERL 2256-ZZL0820 (Union Carbide Corp. ); "B" stage of impregnated fiber: 45 min at 93° C, Mylar cover; 

cure cycle under 35-N/cm^ (50-psi) pressure: 2 hr at 82° C, 3 hr at 148° C. 

^Epoxy compatible (Union Carbide Corp. ). 

Proprietary (Hercules Corp.). 

Poi - Owens Corning Fiberglass Co. 

Composites were fabricated by means of a unidirectional layup of a number of 
B -staged plies to yield the thickness desired. Most of the composites consisted of 
fibers of one particular type. Some hybrid composites were also fabricated that con- 
sisted of two fiber types in the layup with selected thickness and position of each. The 
composites were cured under heat and pressure in a matched-die mold. Complete curing 
conditions are included in table VII. 

Miniature Izod specimens were machined from the fabricated composites in both the 
longitudinal and transverse directions. The finished specimen dimensions were 7. 9 by 
7. 9 by 37. 6 millimeters. 


Test Apparatus and Procedure 

The impact machine used was a modified Bell Telephone Laboratory pendulum type 
(fig. 9). The design capacity of the pendulum was 240 centimeter-newtons (27 in. -lb). 
Addition of weights to the pendulum increased the capacity to 1010 centimeter-newtons 
(114 in. -lb). The striking velocity of the pendulum was 345 centimeters per second. 
The Izod specimens were struck at their free end, 22 millimeters from the edge of the 
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Pendulum release 



Figure 9. - Miniature Izod pendulum-shaped testing machine. 


grip. The specimen length in the grip was 14 millimeters. A "dead weight” load was 
applied to the grip to assure uniform gripping of specimens. 

Composites of one particular fiber were tested in both the longitudinal and trans- 
verse directions. Hybrid composites were generally tested in the longitudinal direction 
with the plies parallel to the striking pendulum. The angular displacement of the pendu- 
lum after impact was an inverse measure of the impact energy. Typical fractured 
specimens from this method of testing are shown in figure 10. 
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S -glass PRD-49 




HTS/Thorne!-50S HMS/Modor-I 



Hybrid composites 


HMS/PRD-49 



071-746 


Figure 10. - Miniature izod impact specimens showing various modes of 
failure. 


EXPERIMENTAL RESULTS AND DISCUSSION 
Longitudinal and Transverse Impact 

Several specimens of each composite system were tested in longitudinal and trans- 
verse impact; also specimens from the matrix system were tested. The results are 
presented in figure 11. The scatter is indicated by the light lines within the bar. 

Photomicrographs of typical fracture surfaces are shown in figure 12. Note the 
fracture modes, cleavage, and cleavage with fiber pullout. (Photographs of the frac- 
tured specimens are shown in fig. 10. ) The observed fracture modes for various com- 
posites are summarized in table VIII. Impact resistance is plotted against short-beam 
intralaminar shear strengths for several of these composites in figure 13. The intra- 
laminar shear strengths are needed to assist with the theoretical impact resistance rank- 
ing of the test specimens. 

Measured results of longitudinal impact normal and parallel to the lamination direc- 
tions were identical. This is to be expected in unidirectional composites with nonmetal- 
lic matrices. 
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Figure 11. - Miniature Izod impact energy of fiber/ERL 2256-ZZL0820 composites. 
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tc) Longitudinal fracture of HTS composite. 


(d) Transverse fracture of HTS composite. 


Figure 12. - Scanning electron micrographs of fracture surfaces of graphite composites resulting from impact load. X600. 
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TABLE VIII. - SUMMARY OF OBSERVED FRACTURE MODES ON 


MINIATURE IZOD IMPACT TEST SPECIMENS 


Fiber; epoxy com- 

Longitudinal 

Transverse 

posites 

Brittle- 

Debonding 

Delamination 

Brittle- 

Fiber 


ness 

plus fiber 


ness 

splitting 



pullout 




HTS 

Yes 

Little 

— 

Yes 

— 

Thornel-50S 

Yes 

Some 

— 

Yes 

Yes 

Thornel-50 

--- 

Yes 

Yes 

Yes 

Yes 

(polyvinyl alcohol) 






Modmor-I 

--- 

Yes 

Yes 

Yes 

--- 

S-glass 

Yes 

Very little 

One layer 

Yes 

--- 

PRD-49 

--- 

--- 


--- 

--- 

HTS Thornel-50S. HTS 

Yes 

Yes 

Yes 

Yes 

--- 

HMS Modmor-I HMS 

Yes 

Yes 

Yes 

--- 

--- 

HMS PRD-49 HMS 

Yes 

Yes 

Yes 

Yes 

— 


1 . 6 ,— 


^ l.2\- 


o 

< 
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INTRALAMINAR SHEAR, PSI 

Figure 13. - Experimental results of transverse impact energy as a func- 
tion of intralaminar shear for various fiber/resin composites. 



Discussion of Experimental Results and Comparison of Ranking 


Examination of table VIII and figure 11 reveal the following: 

(1) Those composites which exhibit more than one fracture mode have higher impact 
resistance in general. 

(2) Composite transverse impact results in brittle fracture and the value is consid- 
erably lower than that of the matrix. Some fiber splitting occurs in the Thornel-fiber 
composites. 

(3) The hybrid composite experienced two or more fracture modes. 

Averaged values of the experimental results are summarized in table VII. The last 
four columns of this table contain the ranking with respect to impact resistance. As can 
be seen, the measured and predicted ranking is identical. The predicted ranking was ob- 
tained as follows: For the longitudinal ranking, equation (13) was used in conjunction 
with table I and figure 13. For the transverse ranking equation (6) was used in conjunc- 
tion with figure 13. The use of figure 13 for the transverse strength is acceptable be- 
cause both intralaminar shear and transverse composite strengths exhibit similar 
trends. 

It is interesting to note in table VI that one of the hybrid composites had larger im- 
pact resistance than either of the two constituent composites. The explanation is that the 
hybrid composite had more delaminated surfaces. This, of course, is the essence of 
the hybrid composite concept for improved impact resistance. 

The important point to keep in mind from this discussion is that theoretical expres- 
sions can be constructed to predict impact resistance at least on a qualitative basis. 
These expressions can be used to guide research for constituent materials and design 
concepts for improved impact resistance. 


CONCLUSIONS 

Results from this investigation of gross-type-impacts of composites involving rela- 
tively long impact contact times lead to the following conclusions: 

1. The impact resistance of unidirectional composites is ranked using elementary 
composite mechanics, and criteria are presented to guide design for improved resist- 
ance. 

2. Theoretical results show that, in composites with high fiber-to-matrix modulus 
ratios, the longitudinal impact resistance is fiber controlled ; When this ratio is twenty, 
the matrix contribution is less than 5 percent. However, the transverse and shear im- 
pact resistances are matrix controlled. 
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3. Theoretical results show that, in composites with fiber-to-matrix modulus or 
strength ratios of about four, the longitudinal impact resistance could be matrix con- 
trolled. In this case, the presence of microresidual stresses decreases the impact 
resistance considerably. 

4. Theoretical considerations indicate that the impact resistance can be improved 
by designing the composite so that fiber breakage, fiber debonding with fiber pullout, and 
partial delamination take place at the same time. Any combinations of these fracture 
modes will also improve the impact resistance. 

5. Theoretical considerations also show that the impact resistance is sensitive to 
void and fiber contents and to certain fabrication factors which are reflected in the in- 
situ constituent properties. 

6. The experimental results indicate three prevalent longitudinal failure modes due 
to impact. These are cleavage, cleavage with some fiber pullout, and cleavage com- 
bined with partial delamination due to intralaminar shear failure. 

7. The transverse failure mode was cleavage. The fracture surface included ma- 
trix fracture, fiber debonding, and some fiber splitting. The experimental results 
showed that the impact resistance was the same whether the specimen was impacted 
parallel or normal to the lamination direction. 

8. Ranking of predicted results was in good agreement with that of measured results 
from notched Charpy impact, cryogenic fracture toughness, stress intensity, and un- 
notched Izod impact. 

9. The hybrid composite concept is an efficient composite design to combine high 
strength and high stiffness with high impact resistance. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 18, 1971, 

129-03. 
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APPENDIX - DETAILED DERIVATIONS OF GOVERNING EQUATIONS 


The detailed derivations which lead to the equations presented and discussed in the 
section THEORETICAL INVESTIGATION are as follows: 

Longitudinal impact (see fig. 1(a)): 


tj _ 1 c e * v = Zl * T V 

U L - g Z11T Z11T 


2E 


Z 11 


From micromechanics (refs. 7 and 8) 


E Z 11 " ' V E f 


E 


k + (l - k f ) -5“' 
E f 


(Al) 


(A2) 


S Z 11T “ ^ " V S fT 


k f%T + /3 m (0 - k f ) — 
E f 


(A3) 


Substituting equations (A2) and (A3) in equation (Al), neglecting terms E m /Ej, and 
simplifying yield 


IEd _ U L _ C 1 k v ) k f / 3 fT S f T 

V 2E f 

Transverse impact (see fig. 1(b)): 

S 2 

tj - 1 c y _ Z22T -y 

U T “ 7 b Z22T e Z22T V " ZZ V 

2 ZE l22 

From micromechanics (ref. 8) 


(A4) 


(A5) 


; _ e mpT F 

’Z 22T “ ^22T ~ H 22 

V/x 22 


(A6) 
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Using equation (A6) in equation (A5) and simplifying result in 


IED 


U T ^ 1 / /3 22T g mpT \ 2 

V 2 \^p22 / 


22 


Shear impact (see fig. 1(c)): 


U c 


= -S 


Zl2S e U2S 


V = 


712S 


2G 


V 


1 12 


From micromechanics (ref. 8) 


o _ a e mpS n 
h 12S “ ^12S ~ W 12 

VVl2 


Using equation (A9) in equation (A8) and simplifying yield 


IED 


U S _ l / 8 12S £ mps \ 2 

V 2 \ ^pl2 / 


G Z 12 


(A7) 


(A8) 


(A9) 


(A10) 


Debonding Contribution 

The work done to pull out broken fibers a distance l cr is given by (refs. 

3 and 4) 


U FPO 



XT xy 7rd f dx 


(All) 


U FPO " g N 


fD’^xy* cr 


(A12) 
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From force equilibrium and assuming uniform shear, we get 


1 cr T xy" u f 


nd f = -i 7rd^S fT 


l 


cr 


1 

4 


S fT 

T xy 


Using equation (A13) in equation (A12) and simplifying yield 


(A13) 


U 


FPO 
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(A14) 


where r xy is the shear interface bond strength. It is generally accepted in the compos- 
ites community that j^g i s a measure of interface bond strength. Letting 
T xy = jgg m equation (A14) results in 


S 2 

U FPO = ^fD nA f ~ 


32S; 


(A15) 


1 12S 


Multiply and divide equation (A15) by the area of the fractured surface as follows: 


U 


FPO 


N fE) 7rdj? S 2 T A 


A 32S 


(A16) 


t 12 S 


By definition 


N f D » d ? 

4A 


= u - V k f 


(A17) 


Using equation (A17) in equation (A16) and dividing through by the area give 


IED = UfPO = (1 - -51- 


Ad, 


8 S; 


(A18) 


1 12S 
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Using the definition for j^g from equation (A7) in equation (A18) results 
sired result 


3 l \ 

8 \^12S e mpS G 112/ 

Delamination Contribution 

The energy expanded in delaminating several interply layers of area Aj 
length Ip is given by 


N 


LD 




3 £ 12S 

h 12 , 


^ A D Z D^i 


Assuming the sum in equation (A20) is independent of i yields 


U D ~ g ^ N LD A D^D^ g 


£ 12S 


£12 


and the impact energy density is 


IED = = N LD 


a d*d 


2G 


£ 12 


Using equation (A9) for ^g in equation (A22) and simplifying results in 


IED = ± N 


LD 


^12S g mps \ 2 

^JLtl2 / 


J l 12 


in the de- 


(A19) 


over a 


(A20) 


(A21) 


(A22) 


(A23) 


34 



Combined Longitudinal and Fiber Pullout 


The energy expended to produce this combination of modes is a combination of 
equations (A4) and (A16) each multiplied by the appropriate volume; that is, 


U = 


~ U ' M ♦ f a - V 

2E f 8 \ /3 12S e mpS/ 


/V% 12^ 

H 

i 

V 12S e mpSy 

G z 12J 


(A24) 


Factoring out (1 - k v )L c A^Sj T /2E^ 


and dividing both sides by L C A yield 


IED = = - (1 - k ) — 

T A O V tt, 


L c A 


E. 


q \ ij c/V12S € mp§/ u l 12J 


(A25) 


Note that in equation (A25) refers to the fiber volume ratio for the whole cross sec- 
tion, while kjp refers to the volume ratio of the pullout fibers. 


Combined Longitudinal Fiber Pullout and Delamination for a 
Rectangular Cantilever (See fig. 8) 

The energy expended in delaminating several layers simultaneously for a cantilever 
is given by 



LD 


,b l 


D S Z12S 


(A26) 


Assuming simultaneous tensile failure and delamination near the neutral plane of the 
cantilever yields 


N LD S £ 12S*D b " \ hbS £ 11T 


and 
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(A27) 


*D = 


bS Z 11T 
4N LD S Z 12 S 


Using equation (A27) in equation (A26) and simplifying yield 


U = 


bh 2 11T 


32N LD S £ 12S 


(A28) 


the total energy expended to fracture the cantilever in combined modes is given by 

U = U FLEXURAL + U SHEAR + U FIBER PULLOUT + U DE LAMINATION ^ A29 ^ 


The flexural and shear energies are given by 


U FLEXURAL + U SHEAR 
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(A30) 


Carrying out the integration in equation (A30), requiring that the cantilever will fail 
first by tension, and simplifying result in 
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(A31) 


The shear contribution (second term in eq. (A31)) will be less than 3 percent if 


M 2 > io fill 

G f 12 


(A3 2) 


For a hybrid composite where the shell and core are made from different compos- 
ites, equation (A31) is approximated by 
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(A3 3) 


U FLEXURAL + U SHEAR “ \ 



where 
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(A3 4) 


and 
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(A3 5) 


The notation in equations (A33) to (A3 5) is as follows: superscript a refers to averaged 
properties; k is a ratio; superscripts and subscripts s and c refer to shell and 
core, respectively. 

The fiber pullout energy is given by equation (A15). Using equations (A33), (A15), 
and (A28) in equation (A29), simplifying, and rearranging result in 
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(A36) 


subject to 
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(A3 7) 
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Equation (A37) supplies the condition to satisfy the assumption of either simultaneous 
tensile failure or simultaneous compressive failure and delamination. Equation (A33) 
can be expressed in terms of constituent material variables by using the micromechan- 
ics definitions for strengths and moduli. 


Longitudinal Impact Resistance When Matrix Is Controlling and 
Effects of Microresidual Stresses 

This case arises when E^/E m < 10. Fiber/metallic matrix composites usually 
meet this requirement. The energy stored in the composite is given by (refer to 
fig. 1(a)) 


tj _ 1 g r * y —l HT 


(A38) 


The micromechanics equation for when the matrix controls the failure is given 

by 


S l 11T " ^ ' V S mT 
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Using equation (A39) in equation (A38) and simplifying yield 


ied = 2 . = 


2E 
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(A3 9) 


(A40) 


The presence of microresidual stresses will affect S m ,p. This effect will equal the 
magnitude of the residual stress. The available matrix strength for resisting impact is 


®mT 


= S. 


mT °mR 


(A41) 


where a „ is the microresidual stress in the matrix. 
mK 

methods, it can be shown that 


Using strength of materials 
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Substituting for S mT in equation (A40) and using equation (A42) result in 
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(A43) 


Corresponding equations for the fiber are obtained in a similar manner. The result will 
be analogous to equation (A43) with subscripts f and m interchanged. 
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